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Abstract. Suppose granular material is shaken vertically with z(t) = 
Aocos(uJot). Can we expect to find convection if Aqui^ < gl By means 
of theoretical analysis and computer simulation we find that there is no 
critical T = \Ao\ujQ/g for the onset of convection. Instead we propose a 
modified criterion which coincides with T = 1 for small frequency ujq. 



1. Introduction 

In vertically shaken granular material one observes various macroscopic ef- 
fects, such as convection [1-3], surface fluidization [4-6], spontaneous heap 
formation [7,8], surface patterns [9-12], oscillons [13] and others. The com- 
mon feature of all these effects is that particles change their position with 
respect to each other. If a single hard particle is shaken, one agrees that it 
cannot move as long as during the entire period of shaking the acceleration 
due to gravity g exceeds the upwards acceleration of the shaker. For the 
case of sinusoidal motion z(t) = Aq cos(u;oi) with amplitude this condition 
is given if V = AqlUq/ ' g < 1. 

We want to discuss the problem whether in a granular material particles 
at the surface can separate from each other even if V < 1. At first glance 
the problem seems to be trivial, however, one has to take into account that 
there is a nonlinear interaction of touching grains. 

There is a controversial discussion in the literature whether there is 
a critical value of Froude number r c = |j4o|wq/<7 below which the above 
mentioned effects vanish. In many experimental observations (e.g. [1,5,8, 
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9,14,15]) and computer simulations (e.g. [15]) such a critical number T c 
was found. Several authors believe that the value is T c = 1. In numeri- 
cal simulations, however, some authors have found surface fluidization and 
convection for T < 1 [3,6,16]. Gallas, Herrmann and Sokolowski [17] con- 
cluded already 1992 from their numerical results ". . . that the current belief 
that T determines the degree of fluidization is incorrect" . 



2. A one dimensional continuum model 



We discuss the response of granular material with respect to vertical oscilla- 
tion in the limit of a one dimensional approach: the lowest bead of a vertical 
column of N identical spherical beads is shaken with zq = Aq cos uj^t and 
the other beads move due to their interaction force and gravity g. Adjacent 
spheres interact with their next neighbors by the force 

F k ,k+i = -v^R (/^feffe+i + ai ktk+1 y/£k,k+i) > (!) 

with p and a being elastic and dissipative material constants, i.e. functions 
of Young modulus, Poisson ratio and dissipation rate (for details see [18]). 
£ is the overlap 2R — \z k — Zfc+i| of adjacent spheres of radius R and mass 
m at positions z k , z k+ \. The height of the column is L = 2NR. The net 
force experienced by the k-th bead reads, therefore 

Fk,k+i — Fk-i,k = m 'Zk = 

€ k -i,k ~ Zk^k+i +^R~ a ik-i,kyj £k-i,k - L,k+iV£k,k+i (2) 

Introducing new coordinates u k = z k — 2Rk (k = . . . N) the overlap of 
two adjacent spheres can be written as 

£k,k+i = 2i? - z k+1 + z k = u k - u k+l = u(2Rk) - u(2R(k + 1)) (3) 
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In (3) the discrete value u k has been replaced by u(2Rk) which allows for 
a Taylor expansion. Eq. (2) turns into 
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A and B can be evaluated by Taylor expansion 
d 
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Eqs. (9) and (11) contain only local variables, therefore we can drop the in- 
dex k and the equation of motion (5) in continuum approximation including 
gravity reads 
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(13) 



Eq. (13) defines the abbreviations k and (3. We are interested in the 
critical parameters of driving {Aq, uq) when the iV-th particle loses contact, 
i.e. when ujy > u/v-i- According to the nonlinear interaction of the particles 
the motion of all other spheres is not sinusoidal anymore, instead one finds 
a superposition of many frequencies. We define the "response" R(ujo) as the 
ratio An/Aq where An is the amplitude of the N-th particle at frequency 
too which can be calculated by convoluting the motion ZN(t) with exp(iuot) 
and Aq is the amplitude of the driving vibration. Suppose ANio N /g > 1 the 
N-th particle separates from the — 1-st. If we would find A$ < An the 
critical Froude number T c = Aquoq/q would be less than 1. We will show 
that there is a range for uoq where this is the case. 



3. Solution of (13) in the limit of frictionless motion 

The boundary condition for Eq. (13) is (du/dz) z= L = 0. To find the strain 
of the material without external forcing (shaking) we consider the limit of 
no damping (a = 0). With 
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Eq. (13) turns into 
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Eq. (15) is denned for x 6 [0, 1], its time independent solution U(x) is 

3 



(15) 



U(x) = | 7 2 (xf - l) 



(16) 



The solution of (15) can be considered as a superposition of the static 
solution (16) and a perturbation w(x,t). Inserting u = U + w in (15) yields 
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We are interested in standing wave solutions w = T(r, Q)X(x, fi) of Eq. (17) 
and obtain 
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with O being a real number. Obviously for T(r, $7) one gets T ~ exp(i^r). 
The solution of the spatial equation 



2dx 

can be found using the Ansatz 
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X(x,Q)=x5/(y), y=-V6Qx§ 



which yields 
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Eq. (21) is the Bessel equation of order 2/5. Hence, the solution of (19) is 
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An expression containing J 2 /5 would be a solution too, however, it does not 
satisfy the boundary condition of (15). The factor has been chosen in order 
to assure X(0, fi) = 1. 

4. The criterion for the onset of surface motion 

The above defined response R is the ratio A^/Aq for given driving frequency 
ujq, or S7o, respectively. Since the zeroth particle corresponds to x = 1 and 
the 7V-th to x = we can write 

The response R is an amplification factor, therefore, the value g/R{££) is 
the critical acceleration of the driving vibration. R is larger than 1 for all 
driving frequencies to. This means that for any driving frequency loq and 
driving amplitude Aq the amplitude of the top particle of the column A^ 
at frequency ujq will be larger than Aq. Hence, for Aj\ru> 2 /g = 1, i.e. when 
the iV-th particle separates from the N — 1-st, we find A^uj^/g = T c < 1. 

Therefore, we replace the condition T > 1 which was supposed to be 
the precondition for surface fluidization, convection etc., by 

^l = T>R- 1 M • (24) 

The function i? _1 (o;) over uj is drawn in Fig. 1 (full line). The curve 
reveals pronounced resonances at Eigenfrequencies ojk where R~ l becomes 
minimal. All experiments on surface fluidization and convection which can 
be found in literature have been performed far below the first resonance. 
Therefore, of particular interest to practical purposes is the limit of small 
frequency uq, i.e. below the first Eigenvalue. The Taylor expansion of -R -1 (f2) 
for small Q yields 



r- 1 = i-^n 2 + o(n 4 ) = i-l(^Y u 2 + o{u A ). (25) 

5 5 \gn J 

Given the container vibrates with frequency ojq. Then for the critical 
amplitude Aq of the vibration when the top particle separates, i.e. when 
the material starts to fluidize one finds 

g_ _ 2 ( \ 1 

ujI 5 \gn 



A * = -2--A—2) ■ (») 



Surprisingly even for very small frequencies where R 1 —> 1 one finds that 
the critical amplitude is reduced by a constant as compared with g/u 2 ,. 
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Figure 1. The reciprocal response RT 1 over driving frequency u>o without damping, 
Eq. (23) (full line). The dashed lines show numerical results based on integration of 
Eq. (2) for damping (3 = 0.01 and (3 = 0.03. 

5. Numerical solution including damping 

Eq. (22) describes the behavior of the column of grains for the case of 
purely elastic contact (a = 0). We have not been able to solve the full 
Eq. (13) including damping in closed form, therefore we have solved Eq. (2) 
numerically for different values of damping. The dashed lines in Fig. 1 
display the reciprocal response over fio with fixed amplitude A = 
0.1 mm, elastic constant k = 2000m 7 / 2 /s 2 and L = 50 m. Fig. 1 shows 
that for small frequency and small damping a the theoretical curve (full 
line) agrees well with numerical data. 

Hence, we believe that at least for small enough damping the results 
of the previous section remain correct, i.e. there is a region where one can 
find surface fluidization, convection and other macroscopic effects even if 
the maximum acceleration lies significantly below g. 

6. Conclusion 

We derived an equation of motion for a column of spheres on a vertically 
vibrating table, representing a granular material in one dimensional approx- 
imation. It is shown that the sphere on top of the column TV can separate 
from the N — 1-st even if the table oscillates with \Aq \ uj^/g = T < 1. In 
agreement with numerical simulations by Gallas et al. [17] we came to the 
result that T is not the correct criterion for classification of phenomena 
which occur in vibrated granular material. 

We could show that instead of the widely accepted condition Aqujq /g > 
1 one has to satisfy A^uj^jg > R^ 1 where is a function of uj which is 
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always less than one. Numerical calculations with low damping agree well 
with analytic results. 

The described result is in contrast with several experimental investiga- 
tions where a critical Froude number T c > 1 has been measured. Whereas 
the Froude number is certainly the proper criterion to predict whether a 
single particle will jump on a vibrating table we suspect that this number 
is not suited to be a criterion for surface fluidization of a column of spheres, 
and the more not for a three dimensional granular material. 
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